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Abstract 

Rclativistic heavy ion physics studies the phenomena that occur when a very large (in units of 
QCD scale Aqcd) amount of energy is deposited into a large (in units of Aqq^) volume, creating 
an extended in space and time domain with an energy density that is large in units of Aqcd- 
This includes the mechanism by which the energy is deposited (likely a transformation of the 
colliding Lorentz-contracted "gluon walls" into the strong longitudinal color fields); approach to 
thermalization; and the static and dynamical properties of the created quark-gluon plasma. Of 
particular interest is the fate of symmetries (e.g. chiral SUhii) x SUr{3), scale, and discrete P 
and CV invariances) in hot and dense QCD matter. At present, the program at RHIC has entered 
a stage where new discoveries are enabled by high precision of the measurements; moreover, an 
array of new capabilities will soon be available due to the numerous and significant upgrades. 
Very importantly, we will soon have access to unprecedented energies of colliding ions at the 
LHC. In addition, future RHIC runs at low energies, FAIR at GSI and NICA at JINR will make 
possible the studies of QCD matter at high baryon density. I will describe the current status of 
theoretical knowledge about hot QCD, and the ways in which it may be expected to improve in 
the near future. 
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1. Introduction 

Relativistic heavy ion physics is placed at the intersection of three major directions of 
research: i) small x, high parton density QCD; ii) non-equilibrium field theory; and iii) 
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phase transitions in strongly interacting matter. Indeed, understanding the evolution of 
a heavy ion collision requires a working theory of initial conditions, of the subsequent 
evolution of the produced partonic system, and of the phase transition(s) to the decon- 
fined phase. This (clearly incomplete and biased) overview is an attempt to capture some 
of the recent changes and developments in the theoretical picture of these phenomena 
which have been triggered by an intense stream of the new data from RHIC. I will also 
discuss the theoretical expectations for the physics of nuclear collisions at the new energy 
frontier that will soon be explored at the LHC. 



2. Quantum Chromo-Dynamics of strong color fields 

Before the advent of RHIC, it was widely believed that at collider energies the total 
multiplicities will become dominated by hard incoherent processes. The very first data 
from RHIC (for the overview of the data from the initial stage of RHIC experiments see 
[Il2l3l4j l showed that the measured charged hadron multiplicities in Au — Au collisions 
appeared much smaller than expected on the basis of incoherent superposition of hard 
processes. 

Given that any inelastic 
rescatterings in the final state 
can only increase the multi- 
plicit}p~| we have an exper- 
imental proof of a high de- 
gree of coherence in multi- 
particle production in nuclear 
collisions at RHIC energies. 
At high energies, all par- 
tons in the nucleus located 
at the same transverse coor- 
dinate contribute coherently 
to the scattering process, and 
the density of partons in the 
transverse plane introduces a 
new relevant dimensionful pa- 
rameter - the saturation mo- 
mentum Qs |5)6j . The den- 
sity of partons in a nucleus 
with mass number A grows 
as Ql ^ A^/'^, so the sat- 
uration momentum becomes 
quite large, Qs ~ 1 GeV for a 
heavy nucleus at RHIC ener- 
gies. One may thus hope that 

the corresponding coupling constant as{Qs) is reasonably small to make weak coupling 
methods work in describing the dense partonic system. In this case the occupation num- 




Fig. 1. Comparison of the prediction for charged hadron multi- 
plicities within the CGC-based KLN model compared to other 
approaches in central PbPb collisions at y/s = 5.5 TeV per 
nucleon pair at the LHC; from I17| . 



For statistical systems, this is due to the second law of thermodynamics 
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ber of gluon field modes with transverse momenta below Qs is ^ I /as{Qs) and large - we 
are dealing with a semi-classical gluon radiation field 'B', or the "color glass condensate" 
(hereafter CGC). The small x evolution thus occurs in this semi-classical background, 
and is described by non- linear GLR [7], JIMWLK [9] and BK |10|llj evolution equations. 

In a frame where both nuclei move fast (e.g. in the lab frame at RHIC and LHC), we 
are thus dealing with the collision of two dense " gluon walls" containing transverse color 
fields. Immediately after the collision, the flux of longitudinal color fields is stretched in 
between the remnants of the nuclei. This has long been the picture within a string model, 
but in the CGC approach the created longitudinal ^12j ("glasma" [13 ) fields have large 
strength E,B ^ Q^. These strong fields then decay into quarks and gluons creating the 
hot matter through a mechanism the dynamics of which is not yet entirely understood. 
However much progress has been made for example by solving the Yang-Mills equations 
to follow the real-time collision dynamics on the lattice (for concise reviews, see |14ll5j ). 
It is useful to have a simple analytical approach encoding the physics of parton saturation 
- an example of such an approach is the KLN model [16^ that evaluates multi-particle 
production in nuclear collisions as a classical radiation off the colliding saturated gluon 
walls. It has been found to successfully describe RHIC data on hadron multiplicities; 
Figjl] shows the prediction for the LHC - again the predicted multiplicities are smaller 
than in most other approaches. 

3. The strongly correlated Quark-Gluon Plasma 
3.1. Strong or weak coupling? 

The quest for the quark-gluon plasma has entered a new era when RHIC has turned 
on, and will intensify at the LHC. We have already learned that the quark-gluon plasma 
at the temperatures accessible experimentally does not at all resemble a quasi-ideal gas 
of quarks and gluons, and the interactions among its constituents are very strong; these 
interactions for example transform a highly viscous dilute substance into what is perhaps 
the most perfect liquid that ever existed. Nevertheless, here I use the term "strongly 
correlated" instead of more familiar "strongly coupled" because at present we do not 
yet know whether the strong dynamical correlations that exist in the system necessarily 
require the coupling constant to be large, or can be reproduced by weak coupling methods 
through an appropriate resummation. 

At the highest temperatures T accessible experimentally, the strong coupling constant 
is about as(27rr) ^ 0.3 which may justify an attempt to use weak coupling methods. 
On the other hand, the corresponding 't Hooft coupling A = g'^N^ = 47rasiVt, ^ 10 
may be large, as argued by the proponents of the strong coupling, large Nc methods. 
Moreover, the very distinction between weakly and strongly coupled plasmas is not well- 
defined: we know that for different observables the corrections (both perturbative and 
non-perturbative) can be very different, the situation already familiar from T — QCD. 
For example, current correlation functions in the quark vector channel — ip^^ijj reach 
the perturbative behavior already at quite small values of '--^ 1 GeV^ whereas in the 
scalar gluon channel there are strong non-perturbative corrections up to ~ 10 GeV^; 
see e.g. [20]. Below we will see that this situation persists also at finite temperature; the 
question "weakly or strongly coupled" thus is ill-defined unless we specify the observable. 
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3.2. First-principle numerical calculations: lattice QCD 

Lattice QCD is a most valuable source of information about the properties of QCD 
at finite temperature. Recent rapid developments in this field (owing both to the better 
methods and the availability of more powerful computers) have brought us very close 
to the understanding of QCD thermodynamics with physical quark masses. The energy 
density and the pressure of the QCD plasma as a function of temperature |19| is shown 
in Fig [2j left. To realize what temperatures are accessible at RHIC and LHC, we first 
evaluate the initial energy density in central AuAu or PbPb collisions using the values of 
saturation momentum Qs describing the multiplicities measured at RHIC and expected 
at the LHC, and then convert this energy density into the corresponding temperature. 
Following this procedure we get 

Trhic - 300 ~ 400 MeV; T^hc - 500 600 MeV. (1) 

One can see that already at the initial temperature achieved at RHIC the energy density 
is well above the critical value corresponding to the transition from hadronic gas (with 
^ degrees of freedom) to the quark-gluon phase (with ~ degrees of freedom) . One 
can see also that even that at the LHC the energy density and pressure do not yet reach 
their values expected for an ideal quark-gluon gas. Does this mean that the dynamics 
explored at RHIC and LHC will be qualitatively the same? 
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Fig. 2. Left: energy density and three times the pressure of QCD matter with almost physical quark 
masses (two light and a heavier strange quarks) as a function of temperature; right: the behavior of 
(e — 3p)/T'* above the deconfinement phase transition. From 1191 . 

To answer this question we need to look at the lattice data in more detail; consider 
for example the thermal expectation value of the trace of the energy-momentum tensor 
{9'-^)/T'^ — {t — 3p)/T^ that has to vanish for an ideal massless quark-gluon gas. The 
lattice result is shown on Fig |2] right; one can see that this quantity exhibits a steep 
temperature dependence in the range of temperatures given by Eq. ([T]) , and is significantly 
different at the initial temperatures reached at RHIC and the LHC. There are at least two 
reasons why the measurement of (9^) is important. First, in field theory the trace of the 
energy-momentum tensor is equal to the divergence of scale (or dilatational) current 
generated by the rescaling of coordinates, d^s^ — f?^. Therefore the decrease of (e— Sp) /T'^ 
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towards higher temperatures signals the approach to an approximately scale-invariant 
behavior; the issue of scale invariance is crucial in deciding whether conformal theories 
(e.g. maximally supersymmetric A/" = 4 Yang-Mills theory) can provide a qualitative 
insight into the properties of quark-gluon plasma. Second, the very character of the 
deviation from conformal behavior can tell us a lot about the properties of the plasma 
close to Tc- The behavior of (e — 3p)/T* has been fitted (see curves on Fig[| right) by 
the following functional dependence 

/e — 3p\ _ 3 4 b c 

where bo is the coefficient of the [3 function, and b and c are the adjustable parameters. 
The first term of course arises from perturbation theory. The term proportional to 1 /T** 
stems from the zero-temperature non-perturbative physics (for example, within the bag 
model it is proportional to the bag constant B , c ^ B); numerically, it appears relatively 
small PJJ. The second term ~ is the most interesting: it can be thought of as a 

"power correction" to the equation of state arising from an operator of dimension twcp^ 
The physics behind this term [21'I22!231 may hold the key to understanding the plasma 
structure close to T^. such a term can arise from a contribution from magnetic monopoles 
or, in other language, from the strings that may still be present in the plasma SitTi.<T< 
(2 3)Tc. The massless excitations of strings [24 with tension a at finite temperature 
contribute the ~ cr/T^ term to the trace anomaly. Such "remnants of confinement" may 
be responsible for the strongly correlated nature of the quark-gluon plasma at moderate 
temperatures. The role of magnetic degrees of freedom in the plasma close to Tc has been 
emphasized e.g. in Refs. |25l25] . 

Another important piece of 
information on the dynami- 
cal degrees of freedom in the 
plasma comes from the lattice 
measurements of the fiuctua- 
tions of the conserved quanti- 
ties: baryon number, electric 
charge and strangeness [TB], 
see Fig. |3] One can see that 
these fluctuations reach the 
values expected for an ideal 
quark gas already at T ~ 
1.5 Tci in sharp contrast to the 
case of We are thus deal- 
ing with a very different be- 
havior in the vector channel 
(conserved quantities are the 
temporal components of the 



^ Note however that within the current precision this term cannot be unambiguously separated from a 
logarithmic dependence resulting from the running coupling constant. 
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Fig. 3. Quadratic fluctuations of baryon number, electric 
charge and strangeness, normalized to the corresponding free 
quark gas values; from |18j . 



vector current, ~ ipjQ^) and in the scalar one, similarly to the case of zero tempera- 
ture. 

3.3. Perturbation theory and quasi-particles 

At weak coupling g, quarks and gluons are the dynamical degrees of freedom in the 
quark-gluon plasma. The widths T of quark and gluon quasi-particles are small compared 
to the typical thermal scale given by the temperature T, F <C T because they are sup- 
pressed by the powers of the coupling constant g I. Gluons acquire a dynamical Debye 
mass mo ^ gT, and the color Coulomb potential is screened at distances ^ I/tud- Pow- 
erful resummation techniques have been developed to describe the thermodynamics and 
transport properties in this regime; a particularly useful Hard Thermal Loop method |27j 
utilizes the hierarchy of magnetic g^T, electric gT and thermal T scales: g^T <^ gT <^ T. 
Resummation methods allow to extend the applicability of the weak coupling approach 
down to moderate temperatures, perhaps as low as T ~ (2 3)Tc (where Tc is the de- 
confinement transition temperature); for a review, see [28j . The parton scattering cross 
sections in this regime are small, ^ 5^, and this leads to large shear viscosity r] and small 
bulk viscosity ^; dimensionless ratios of these quantities to the entropy density are (up 
to logarithms) |29l30j 

^^^»1; (3) 

Things become much more complicated at temperatures Tc < T < (2 3)rc- In this 
temperature range (most relevant for RHIC experiments) the coupling constant g^{2TrT) 
is not small, and the quasi-particle description breaks down. In other words, the gluon 
clouds of the quark and gluon quasi-particles become so dense and the interactions be- 
tween the clouds of different partons becomes so intense that the representation of this 
system as a superposition of individual dressed partons fails, even as a rough approxi- 
mation to reality. 

3.4. AdS/CFT correspondence, unparticles, and the "unplasma" 

Much of our physics intuition is based on the quasi-particle description - most of 
us have a difficulty to imagine a composite object without well-defined constituents. 
However this kind of behavior can be encountered for example in the studies of conformal 
theories. Conformal theories do not possess dimensionful scales, so the spectral density 
of an excitation should look like a structureless mass distribution - instead of a particle 
with a well-defined mass we are dealing with a broad amorphous " unparticle" |31j . This 
situation is somewhat reminiscent of what we expect to happen to quark-gluon plasma 
at strong coupling - no well-defined quasiparticles exist, just very broad "unquarks" and 
"ungluons" . We are thus dealing with what may be called an "unplasma" of such unquark 
and ungluon states. 

While we still do not know how to treat QCD analytically at strong coupling, there 
has been a breakthrough in understanding the dynamics of maximally supersymmet- 
ric conformal J\f = 4 Yang-Mills (YM) theory based on the holographic correspondence 
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between the field theory describing behavior on AD Minkowski boundary and the super- 
gravity in the 5D bulk Anti-de Sitter space, supplemented by a 5D sphere, AdS^ x 
[32133134] . The AdS^ metric is unique if we want to extend Minkowski boundary space 
into the fifth dimension preserving the conformal invariance. A particularly useful feature 
of this correspondence is the duality between the strongly-coupled Yang- Mills theory and 
the weakly coupled gravity that can be treated by semi-classical methods. On the other 
hand, weakly coupled YM theory that is amenable to perturbation theory analysis is 
dual to strongly coupled quantum gravity. In particular, there is an intriguing relation 
between YM dynamics at weak coupling and the dynamics of gravitational collapse at 
strong coupling [35] . Most of the applications however deal with the strongly coupled YM 
theory using the duality to translate some hopelessly difficult field-theoretical problems 
into treatable exercises in classical gravity. The finite temperature T is introduced on the 
gravity side by putting a black hole in the center of the AdS^ - note that the Anti-de 
Sitter metric is the solution of Einstein equations for a negative cosmological constant, 
and so " gravitates onto itself . The temperature T and the entropy density s of the YM 
plasma are thus associated with a Hawking temperature of the black hole and the entropy 
associated with its event horizon. 



3.5. Shear and bulk viscosities; plasma as a "perfect liquid" 

As a specific example, let us mention an elegant calculation of shear viscosity carried 
out in Ref. (36j . The shear viscosity is defined by Kubo's formula through the low- 
frequency behavior of the correlation function of the energy-momentum tensor. On the 
Minkowski boundary, the energy-momenum tensor is an operator composed of YM fields. 
In gravity, the energy-momentum tensor excites gravitons, so the problem of computing 
the correlation function reduces to evaluating the amplitude of graviton propagation in 
the AdS^ background. We are interested in solving this problem in the regime of strong 
YM coupling, and thus of weak gravitational coupling. In this case the problem reduces to 
evaluating the absorption cross section aabs of the graviton by the black hole, determined 
by the area of the horizon; since the entropy of the black hole is also given by the area 
of the horizon, one arrives at a dimensionless bound rj/s = I /At: that is achieved at 
strong coupling. There is an ongoing debate on whether or not this bound is universal. 
Addressing bulk viscosity requires introducing deviations from the conformal behavior; 
however this problem may still be treatable using the gauge-gravity correspondence (see 
e.g. [37] and references therein). 

In heavy ion collisions, an economical way of describing the evolution of the system 
is provided by relativistic hydrodynamics, which transforms the gradients of the initial 
parton density into the momentum flow of the produced hadrons. Of particular interest 
is the azimuthal asymmetry of the hadron momentum distribution parameterized by the 
"elliptic flow" V2 (for a review, see [38 ). This quantity appears sensitive to the value 
of shear viscosity, as illustrated in Fig. |4] left, [33] ■ It appears that the data favor small 
values of shear viscosity that are not far from the "perfect liquid" bound discussed above. 
Viscous effects contribute not only at the quark-gluon plasma stage of the evolution, but 
also after the hadronization; these effects can be introduced through a hadron cascade, 
see Fig. |4] right [40]. 
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4. Topology-induced parity violation 
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Non-Abelian nature of QCD 
has a very important impli- 
cation for the structure of 
the vacuum state. Indeed, 
the classical Yang- Mills equa- 
tions (non-Abelian analogs of 
Maxwell equations) possess 
non-trivial vacuum solutions 
- instantons - corresponding 
to the mapping of SU{2) sub- 
group of the gauge symme- 
try SU{3) onto the group 
of three-dimensional rotations 
S3 [H]. Instantons thus cou- 
ple rotations in space to ro- 
tations in the space of color. 
In Euclidean space-time, in- 
stantons are static localized 
objects. In Minkowski space- 
time, instantons describe tun- 
neling transitions between the 
states with different topolog- 
ical Chern-Simons numbers 

[12] v of the SU{2) ^ S3 mapping |43I44I45] . At finite temperature, the transition be- 
tween the vacuum states with different topological numbers can occur not only through 
quantum tunneling, but can also be induced by a classical thermal activation process, 
through a "sphaleron" [46| that induces baryon number violation in electroweak theory 



10 20 30 40 50 60 70 80 



Fig. 5. Correlations of charged particles with respect to the 
reaction plane that reveal parity violation in AuAu and CuCu 
collisions at y/s = 200 GeV per nucleon pair; from 1521 . 
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and helicity non-conservation in QCD |47| . The rate of the sphaleron transitions T^ph is 
not exponentially suppressed, and is proportional to T^ph ~ (with a numerically 

large coefficient [48^). Sphalerons describe a random walk in the topological number; in 
a volume V and time period T we get the topological number (i/^) = TsphVT. 

In off-central collisions, heavy ions possess a very large relative angular momentum 
and create very strong magnetic fields. In this situation, the presence of topological 
charge was predicted [IS] to induce the charge separation with respect to the reaction 
plane (and thus the electric dipole moment of the produced quark-gluon plasma). Soon 
afterwards, an experimental observable sensitive to the effect has been proposed [50] . 
The first preliminary results have been reported by STAR Collaboration in Ref. [51] : 
Figj5]shows the reported in [52] preliminary results for the charge asymmetries measured 
in AuAu and CuCu collisions at RHIC. A detailed theory of the interplay between the 
chiral charge and the background magnetic field responsible for the charge separation 
(the "chiral magnetic effect") has been developed in Refs. [53I54|55] . The generation 
of chirality in the quark-gluon plasma that is responsible for the charge separation is 
the QCD counterpart of the generation of baryon asymmetry in the electroweak phase 
transition in the Early Universe [56157] . 

5. Summary 

The theory of hot and dense QCD matter is developing very fast, stimulated by the 
stream of high quality data from RHIC and the anticipations for the LHC. It is very hard 
to review all of the developments within a limited space; the most glaring omission in 
this talk is perhaps the theory of hard probes interacting with the QCD medium. I have 
no doubt that the studies of jets, high transverse momentum particles, quarkonia, heavy 
quarks and electromagnetic radiation that proved so important at RHIC will provide a 
lot more information about QCD matter at LHC. Another omission is the physics of 
QCD matter at high baryon density that will be explored at low-energy runs at RHIC, 
FAIR and NICA - a topic of fundamental interest. To summarize, we have a lot to look 
forward to - this is an exciting time to be a heavy ion physicist. 

I thank Itzhak Tserruya and his colleagues for their invitation to Eilat and the impec- 
cable organization of the PANIC '08 Conference. 
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